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LONG-RANGE ATOMIC INTERACTION POTENTIALS FOR NITROGEN 


By Willard E. Meador 
Langley Research Center 

SUMMARY 

Direct correlations between molecular and atomic wave functions are employed to 
find the coefficients of the two-electron exchange energies in the valence-bond descrip- 
tion of the interatomic forces between colliding nitrogen atoms. Whereas the X *2* 

and '2* potentials are identical with those deduced previously, the A ^2* and ^2* 
potentials are very different. Not only are the long-range tails of the second pair 
revealed as mirror images of each other, as are the first two, but also each is expressed 
as a simple fraction of the lowest curve. A more general discussion of the different 
types of diatomic bonds and antibonds, including comparisons with those of the current 
literature, also is presented. 


INTRODUCTION 

The long-range tails of interatomic potentials are important for many calculations 
of the nonequilibrium properties of collision-dominated gases. Often used for this pur- 
pose is the approximation of perfect pairing (ref. 1), hereafter called APP, in which 
molecular energies are expressed as linear combinations of bond and antibond energies 
based upon symmetry and spin pairing considerations. Although the absolute results for 
individual collisions may be inaccurate, the method is thought to establish fairly reliable 
and useful relations between the various potentials operating in a gas. (See ref. 2 for a 
discussion.) Hence, if a small number of such potentials can be computed accurately 
(see refs. 3 to 5 for ab initio quantum mechanical calculations of the N 2 ground state), 
the others will follow from simple APP theory. 

A principal objective of the present research is to supply the APP bond and anti- 
bond energies for collisions between two ground-state nitrogen atoms and corresponding 
to molecular electron configurations determined from correlation studies similar to that 
of Mulliken (ref. 6) for the N2^2* j state. More general APP multielectron energies 
also are considered and several important relations are discussed. The correlation 
analysis forces the molecule to separate into atoms with the desired spin and angular 
momentum quantum numbers, a requirement which has not always been satisfied in pre- 
vious calculations. Vanderslice, Mason, and Lippincott (ref. 7), for example, assume 



that of the several different electron configurations which correspond to a given set of 
molecular quantum numbers, only the strongest attractive or the least repulsive one dis- 
sociates into normal atoms. Such assumptions are seldom valid and frequently yield 
incorrect results for the interaction potentials. 

The long-range potential tails thus computed for the N 2 ^X N 2 ^A 

N 2 (^Sg), and N 2 ^ 2 ^ states are compared with each other and with Gilmore’s extrap- 
olations (ref. 8) of Rydberg-Klein-Rees (designated "RKR”) data. Although the latter 
curves are somewhat uncertain in the region of interest, a pronounced preference is 
shown for the present calculations over those of reference 7 mentioned previously. In 
addition, the present theory yields much simpler relations between the energies than does 
reference 7 because each of the three excited potentials can 'be expressed for the first 
time directly in terms of the molecular ground-state energy. Vanderslice, Mason, and 
Lippincott, on the other hand, are left with two independent potentials, one of which they 
obtain by a delta-function computation of N2^2*^ and the other by a rather question- 
able analysis of predissociation data. The latter procedure is unnecessary in the pres- 
ent refinement. 


SYMBOLS 


a,b 


c b c 2> c 3 


E 0 

E(He) 

F 

H 

I ,i 


atoms a and b; reference to nuclei when used as subscripts; atomic 
functions 

variational coefficients in equation (18) 
absolute electron charge 
ground-state energy of atomic hydrogen 
ground-state energy of atomic helium 
arbitrary operator in equation (Cl) 

Hamiltonian operator 

atomic functions, p x + ipy and p x - ipy 
integers 


2 



J i 

L 


Lmax 

L + ,L_ 

m l 

M s 

O 

°L 

Px’^y’^z 

Qi 


r b 

r 12 

R 

S 


Ti 

x,y,z 


exchange energy for an i-electron bond; subscripts a or it refer to i = 2 

quantum number for total orbital angular momentum 

maximum value of L 

step-up, step-down or ladder operators 

magnetic quantum number for total orbital angular momentum 

magnetic quantum number for total spin 

atomic function p z 

projection operator 

p-type atomic orbitals 

Coulomb energy for an i-electron bond 

ratio defined by equation (31) 

distance between electron i and nucleus a 

distance between an electron and nucleus b 

distance between electrons 1 and 2 

internuclear separation 

overlap integral 

s-type atomic orbital 

kinetic energy operator for electron i 

abbreviated notation for the atomic functions p x , Py, and p z ; also 
Cartesian coordinates 
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spin functions 


«,/ 3 


variational parameter 






a, a 


T i 


molecular orbitals 

7 r-type molecular orbitals 
a-type molecular orbitals 
volume element for electron i is dTi 


<t> 


$ d 




interaction potential; sometimes used with subscripts a (attractive) and 
r (repulsive) 

London dispersion energy 

wave function 


I//}, ^2 5*^3 special wave functions defined by equation (11) 

special wave functions i //j and 2~ ^/ / 2 + ^ 3^5 respectively 

Special notations: 


l v + 3 v + 5 V + 

/j . Zj m Zj « 

g’ u’ g» 




7 s;,a 3 s;,x 


\ molecular states 


j 


4 2 2 

S, D, P atomic states 
< > Dirac integral defined by equation (Cl) 

(xy|z) atomic state function defined by equation (B2) 

1,2, 3, 4 sets of spatial coordinates for electrons 1, 2, 3, 4 

Primes are used as indices to distinguish different functions. 
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ATOMIC AND MOLECULAR CORRELATIONS 


With the z-axis defined as coincident with the line joining the two nuclei of a nitro- 
gen molecule, the conventional molecular orbital (MO) notation of a, a*, ir x , rr*, 7 r y , 
and 77* is chosen for the individual states accessible to the electrons of a p3 configu- 
ration on each atom. The various two-electron bonding and antibonding possibilities are 
reviewed and summarized in appendix A, it being assumed that only the interactions 
between electrons in similarly oriented p-type atomic orbitals (AO) will contribute to the 
interatomic potentials. Overlap between the atomic orbitals of different atoms is the 
essential consideration. 

As is explained also in appendix A, doubly occupied molecular orbitals such as (o) 2 
must not be regarded too literally at large distances because of the configuration inter- 
action with (a*) 2 ; hence, the proper wave function for (o) 2 is the valence bond (VB) 
expression z a (l) z b (2) + z b (l) z a (2) constructed from p z orbitals z a and z b cen- 
tered, respectively, on the nuclei a and b. Likewise, with the understanding that all 
dissociations of present interest leave one electron in each of z a and z b , the triplet 
relation of equation (A4) is employed always in place of the singlet relation of equation (A5) 
for the notation (cr)(cr*). Similar statements apply to orbitals with other symmetries. 
Although the dissociation products of a given molecular state can be presented in a num- 
ber of ways, corresponding to the various orientations of the individual atomic spins, only 
one set of atomic functions is ever required for the desired correlations. A complete 
listing of the pertinent atomic states is given in appendix B. 

The following dissociations of molecular nitrogen are significant: 


(c) ( 7 T X ) (7T y ) 

— - ( x a y a z a) + (| x b y b z b) 

(i) 

w 2 (*x) 2 (» y)(*$) 

— - ( x a y a z a) + (y b i z b x b) 

(2) 

( a ) 2 ( 77 x)( 7 r x)( 7 I ’y ) 2 

— — ( x a y a z a) + ( x b|yb z b) 

(3) 

(a)(a*)(7r x ) 2 (7T y ) 2 

— ( x a ya z a) + ( z b| x b y b ) 

(4) 

(a) 2 (7r x )(7r*)(7r y )(7r*) — 

— ( x a ya z a) + ( x bybl z b) 

(5) 

( cr) (cr* ) (77 X ) 2 (7T y ) ^ 77 *) _ 

— — ( x aya z a) + (y b z b| x b) 

(6) 

(a)(a*)(w x )(7r*)(7T y ) 2 — 

— ^ ( x a y a z a) + ( z b x b|yb) 

(7) 
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( 8 ) 


7 S+: (a)(a*)(77 x )(7r*)(7r y )(7r*) ^ (x a y a z a ) + ( x b Y b z b) 

E a ch of the parentheses to the right of the arrows is symbolic of determinants defined as 
in equation (B2), the interconnecting positive signs between them referring schematically 
to N2 - N + N. The dissociations (2), (3), (6), and (7) are not deduced from pure molec- 
ular states but rather from mixtures of 2 and A symmetries; the correct combina- 
tions to establish definite 2 states occur automatically in the subsequent analysis. 

It is noted that only the first atomic parentheses in each of the relations (2) to (7) 
represents an observed atomic state, namely, the n[ 4 s(m s = 3/2) of equation (B6), 
whereas the second is only a part of equation (B4), (B5), (B7), or (B8). Accordingly, the 
following linear combinations are required: 

p(x X 2g - 4 S + 4 s) = I// (o) 2 (%) 2 ( 7Ty) 2 (9) 

3 2+ ~ 4 S + 4 s) = 3 -1/ V ( a ) 2 ( 7r x) 2 ( 7r y) (^y) + 3 " 1/2 ^ (° r ) 2 ( 7r x)( 7r x)( 7T y f 

+ 3- 1/2 ^[(a)(a*)(77 x ) 2 ^ 2 ‘ (10) 

^^2g - 4 S + 4 s) = 3" 1/2 ^ W 2 (*x)(*x)(*y)($j + 3" 1/2 ^|(a)(a*)(7r x ) 2 (7r y )^ 

+ 3 “ J / 2 1// J( a) ( a* ) (w x ) ^7r x ) ( tt y ) 2 ® + i// 2 + ^3) ( n ) 

and 

^ 7 2y - 4 S + 4 s) = p ( <J )( cr *)(^x)^x)( 7r y)( 7r yj ( 12 ) 

where, the notations ?^ 2 , and if / 3 are defined in equation (11) for future reference 

and are considered to be orthonormal. 

The coefficients on the right-hand sides of equations (10) and (11) are of equal mag- 
nitude because of the similar situations in equations (B7) and (B8); by way of contrast, 

^ 5 2+ ~ 4 S + 2 d) = 2(6)" 1/2 ^ (*) 2 (ff x )(4)(*y)( ff y) “ 6' 1/2 ^(a)(a*)(7r x j2^ y j^*J 

- 6 - 1/2 V ( o-) ( a* ) ( tt x ) (tt*) ( 7r y ) 2 = 6" 1/2 (2^i - <// 2 " ^3) (13) 

from equations (B4) and (5) to (7). 
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If the further definitions i//j = and 4 / jj = 2 “ ^^^2 + ^ 3 ) are substituted into 
equations (11) and (13), there results 

- 4 S + 4 s) = 3" 1/2 ^ + (2/3) 1/2 i// n (14) 

and 

j//^ 5 S+ - 4 S + 2 d) = (2/3) 1/2 ^ - 3“ 1/2 i// n (15) 

These expressions are identical with those of Mulliken (ref. 6) mentioned previously. A 
disadvantage of this notation with regard to intermolecular forces is discussed in the next 
section. 

Finally, since (o) 2 forms a stronger two-electron bond than either of ^tt x ) 2 or 
^ 7 Tyj 2 because of the greater overlap involved, an elementary consideration of the non- 
crossing rule for potential energy curves of the same symmetry yields the correlation 
diagrams 


5 Sgt + (2/3) 1 / 2 ^: 4 S + 4 S (16) 

and 

5 Sj: (2/3) 1/2 ^ x - 3" 1/2 ^n: 4 S + 2 D (17) 

as the internuclear separation R increases from left to right. A similar procedure can 
be applied to the dissociations of N 2 ^A 

INTERACTION POTENTIALS 

The combinations of electron configurations in equations (9) to (12), together with 
the two-electron bond and antibond energies (assumed to be additive) of appendix C, are 
sufficient to determine the long-range potentials of two colliding n( 4 s) atoms. Equa- 
tion (11) is preferred over equation (14) because the resonance implied by i//jj between 
a single ir x bond and a single Hy bond is not so directly accessible to the simplest 
APP techniques as are the individual functions \j / 2 and i//g. In the final analysis, of 
course, the results must be equivalent. 

One possible source for confusion remains. Since equations (11) and (13) are spe- 
cial examples of the more general variational function 
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= Ci\f/i + C2^2 + c 3^3 


(18) 


and since the variational process decrees that two of the coefficients Cj, C 2 , and C 3 
must vanish unless the energy integrals <^|h|^j> differ from zero when i ^ j, inter- 
actions between the individual configurations of 4 / 2 , and 4 / s cannot be ignored in 

considerations of the total energy. At the large internuclear separations R of present 
interest, however, the sole function of such cross terms is the separation of the energy 
asymptotes corresponding to different dissociation products; accordingly, the interatomic 
forces are not affected. For comparatively small values of R, on the other hand, the 
cross term <4'i\Ji\4 / 2 + ^ 3 > has a negligible splitting effect because of the already large 
energy separation between the states i//j and 4 / u't this statement, in fact, is the explana- 
tion for i//j and 4 ^u existing alone at the far left-hand sides of the correlation dia- 
grams (16) and (17). 

Hence, for the present application, the correct procedure is as follows: first, the 
APP interaction potentials are deduced independently for each substate in equations (9) 
to ( 12 ); secondly, these substate energies are given their proper statistical weights 
according to the squares of the coefficients of the corresponding wave functions; finally, 
the energies are added to yield the total interaction potentials. The results are: 



where J a and J n are exchange energies of the type J 2 in equation (C9). Improved 
estimates would include the London dispersion energy (Van der Waals attraction) in addi- 
tive form (ref. 2 ). 

For purposes of comparison, the energy 

<p(^Z+ - 4 S + 2 I^ = 3 _ 1 (j CT - 4J n J (23) 
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is derived from equation (13). The interaction potential and the reference asymptote both 
differ from equation (21), even though the molecular state is in each case. 

Although the physical origin of the exchange forces is made sufficiently clear by the 
Hellmann- Feynman electrostatic theorem (the enhancement of the negative charge distri- 
bution between nuclei in a bonding situation being related to the exchange forces, which, 
in turn, result from the particular mathematics of the VB procedure for satisfying Pauli's 
exclusion principle), detailed calculations involving equation (C9) are less useful than 
might be supposed. Computations based upon simple isolated hydrogen-like functions are 
too inaccurate. Nevertheless, as was shown first by Vander slice, Mason, and Lippincott 
(ref. 7), much information can be obtained from expressions like equations (19) to (22) 
when they are employed judiciously with experimental data. 

RELATIONS BETWEEN POTENTIALS 


The most obvious characteristic of equations (19) to (22) is that the dependence upon 
J a and is identical in the four states; accordingly, 


4 >(a 3 S+) = -4>( 5 2+) = -3 = 3 _1 4>(x 


(24) 


If the lowest potential curve is provided from ab initio quantum calculations (refs. 3 to 5), 

extrapolations of RKR results, or some other type of experimental data, the others are 

immediately available. Alternatively, delta function calculations (ref. 9) on the repulsive 
7 + 

S u configuration could be used. 

Vanderslice, Mason, and Lippincott (ref. 7), on the other hand, employed only the 
first substate in each of equations (10) and (11) to obtain 



and 



They then deduced the final relations 

4>( 7 £u) = -*(x ^g) 


(25) 

(26) 

(27) 

(28) 
(29) 
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w 


and 



each of which is included in equation (24). Aside from the greater complexity of equa- 
tions (29) and (30), which presuppose the knowledge of two curves for a complete deter- 
mination of the interactions between normal atoms, perhaps the most significant feature 
of the two sets of potentials is the additional symmetry provided by equations (20) and (21). 

Just as X and are mirror images of each other across the dissociation 

& 3 t 5 + 

asymptote, so are A 2 U and 2g. 

As previously explained by means of the correlation diagram (16) and a similar one 
for the A 3 s£ state, a gradual transition must occur from equations (20) and (21) to 
equations (26) and (27) as R decreases; the present formulation, however, is expected to 
provide a better description over most of the range for which the valence bond theory is 
valid since it was designed for that purpose. A simple test of this prediction is furnished 
by numerical computations of the ratio of <f>(x to 4>(a according to three 

different methods: Vanderslice's formulation, the present formulation, and Gilmore's 
RKR extrapolations (ref. 8). 

Because of the similarity of forms between J a and and the fact that each is 
dependent dominantly upon the square of the corresponding overlap integral, can 

be replaced for present purposes with ( of"- Accordingly, the approach of refer- 
ence 7 yields 



for the desired ratio, whereas equations (19) and (20) give the constant value 3. Table I 
summarizes the various results and clearly supports the present technique over a wide 
range of the potential tails if Gilmore's curves are accurate. Vanderslice's formula does 
not agree with Gilmore's potentials where the VB theory is valid. 

TABLE L- RATIOS r OF (p(x TO $(a 3 S+) AS FUNCTIONS 

OF THE INTERNUCLEAR SEPARATION R 


R, 

angstroms 

Ratios r from — 


Equation (31) 

Present paper 

Reference 8 

2.8 

1.04 

3.0 

2.8 

2.6 

1.05 

3.0 

2.7 

2.0 

1.1 

3.0 

1.9 
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5 - 4 - 

An interesting observation concerns the Eg dissociation: Whereas equation (21) 
predicts a purely repulsive interaction, equation (27) states that the potential is attractive 
at large R because of the dominance of J CT over J^. A more realistic representa- 
tion would retain the repulsion in the far regions, depict a maximum in the potential curve 
at a somewhat smaller separation, and culminate in a minimum between 1.5 and 1.6 ang- 
stroms (ref. 8) as the left-hand side of the correlation diagram (16) is approached. It 
appears doubtful that Gilmore's second minimum between 2.8 and 2.9 angstroms is real- 
istic; in addition, he seems to underestimate the Eg curve between 2.2 and 2.9 ang- 

3 + 

stroms in violation of the requirement that the A S u potential should lie midway 
between those of X *2g and 5 Sg. Equations (29) and (30) are invariant in regard to 
the present and Vanderslice theories, and therefore are especially reliable. 


Finally, if the London dispersion energy is added to each of equations (19) 
to (22) and then eliminated between combinations of the resulting expressions, the following 
relations are obtained between interaction potentials which include the effects of 0^: 


and 



(32) 


(33) 


Equation (32) is common also to the analysis of reference 7 with or without 0^, but equa- 
tion (33) is unique to the current model. 


CONCLUDING REMARKS 


Valence-bond expressions have been derived for the outer regions of the four inter- 
action potentials governing collisions between ground-state nitrogen atoms. Although the 
individual potential tails are probably not very accurate in an absolute sense, the rela- 
tions obtained between the four energy states are expected to be reliable. Hence, the 
combination of the present results with a minimum amount of experimental data should 
prove adequate for most thermodynamic and kinetic-theory applications. 

A principal advantage of the present method over the previous and similar work of 
Vanderslice, Mason, and Lippincott is the use of detailed correlations between the molec- 
ular states and the actual atomic products. The results of these correlations are dis- 
tinguished by additional symmetries of the potential tails. Whereas Vanderslice, Mason, 
and Lippincott found only the most repulsive and most attractive potentials to be mirror 
images of each other, the present research predicts the same property also for the inter- 
mediate pair of curves. Finally, the intermediate potentials are expressed for the first 
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time as simple fractions of the ground-state molecular energy and thus require less 
experimental data for the complete semiempirical specification of nitrogen collisions. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., April 15, 1970. 
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APPENDIX A 


THE HYDROGEN MOLECULE 

A brief review of the essential features of two-electron molecular configurations 
(see ref. 10, pp : 340-358, for a more detailed description), especially in regard to the 
differences between molecular orbital (MO) and valence bond (VB) models, seems appro- 
priate for this paper. The interaction energies of H 2 and several other simple mole- 
cules are derived in appendix C. 

With £ = a + b and £* = a - b as the basic molecular orbits, where a and b 
are 1 s-atomic functions (AO) centered, respectively, on nuclei a and b, there are 
possible six different occupations by two electrons: both in £, both in £*, one in | 
and the other in £,* with antiparallel spins, and three orientations of the parallel spin 
combination of one in £ and one in £*. Only four states of H 2 exist, however, one 
for the singlet X configuration and three for the triplet 2 Su, which dissociate into 

normal atoms; accordingly, two of the six MO occupations must yield excited products. 

A study of the configuration interaction between the first two occupations is essential for 
this distinction. 

The expansion of the two-electron MO wave functions £(1) £(2) = U) 2 and 
i*(l) l*(2) = (£*) 2 gives 

U ) 2 + XU*) 2 = (1 - X)[a(l) b(2) + b(l) a(2)j + (1 + X)[a(l) a(2) + b(l) b(2)] (Al) 

where X is a variational parameter and 1 and 2 refer to the electrons. At large inter- 
nuclear separations (negligible overlap between a and b), the variational procedure 
selects X = -1 and X = 1 corresponding, respectively, to the H + H VB function 
a(l) b(2) + b(l) a(2) and the H + + H” ionic state a(l) a(2) + b( 1) b(2). In addition, the 
application of the noncrossing rule reveals the correlation diagrams 

H 2 ^X 1 sjj : U) 2 — valence bond: H + H 

and 

H 2 ^ SgJ: U*) 2 ►ionic: H+ + H‘ 

for the transition from small to large distances. 

The remaining occupations are, with their wave functions, 

H 2f 3s u) : (DU*) ► valence bond £a(l) b(2) - b(l) a(2)J: H + H (A4) 


(A2) 

(A3) 
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and 

H 2 ( 1 sj): (m*) — ionic [a(l) a(2) - b(l) b(2)]: H+ + H~ (A5) 

Configuration interaction is not appropriate for either of the latter states; more specifi- 
cally, the MO and VB descriptions are equivalent at all separations in the relation (A4). 
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APPENDIX B 


WAVE FUNCTIONS FOR ATOMIC NITROGEN 


The atomic eigenfunctions for the total angular momentum (quantum number L) of 
any complex species are obtained most directly by the technique of projection operators. 
With an arbitrary electron configuration in the linear vector space restricted to the total 
magnetic quantum number Ml equal to L, the subsequent operation (ref. 11) by 


L m ax~ L 

0 L = (2L + i): >> 

t=0 


t!(2L + t + 1): 


(Bl) 


will select just that component of the Hilbert space representation of the initial wave func- 
tion corresponding to the desired value of L. Other symbols are L + and L_, the usual 
ladder operators of the electron system, and the maximum value L max permitted for 
L. Operations by L_ are used for Ml < L. 

As an example of the general procedure, the operator O^ is applied to the three- 
electron function 


(IOlO) = (3 1) -1 / 2 


10(1) 

Oo(1) 

0/3(1) 

Io<2) 

Oa(2) 

0/3(2) 

Io<3) 

Oq<3) 

0/3(3) 


(B2) 


where a and /3 are spin states and I = x + iy, I = x - iy, and O = z, the x,y,z 
referring to p-type atomic functions of different orientations. With L max = 2 because 
of Pauli's exclusion principle, the result of the projection operation is 


i// 2 p(m l = l) = 2 -1 / 2 (I0|0) + (il|l) 


= 10 1//2 [2(xy|y) - (xz|z) - 2i(xy|x) - i(yz|z)] 


(B3) 


The remaining functions of present interest for nitrogen are 


* 




d(m l = 0, Mg = 6 _1//2 [2(xy|z) - (zx|y) - (yz|x)] 


2d(m l = 0, M s = -1 


= 6" 1//2 [2(z|xy) - (y|zx) - (x|yz)] 


(B4) 

(B5) 
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and 


i' 


* 

* 


Vm s =| 


4 s(m s = -| 



= 3" 1//2 [(xy|z) + (zx|y) + (yz|x)] 
= 3 _1 / 2 J(x|yz) + (z|xy) + (ylzx)] 



(B6) 

(B7) 

(B8) 


(B9) 


where M s is the magnetic quantum number of the total spin. Step-down operators in 
spin space were employed for the derivation of equation (B5) from equation (B4) and 
equations (B7) to (B9) from equation (B6). 

It is noted that the wave function for n(^p) differs from those of n(^d) and 
n( 4 s) in two important ways: in the first place, the former is complex; secondly, it 
shows double occupations of the various p-atomic orbitals. Accordingly, the corre- 
sponding contributions to the interaction forces between two nitrogen atoms are not those 
of appendix C, but rather must be derived in similar fashion from the ionic functions in 
equations (A3) and (A5) of appendix A. 


* 


% 


16 



APPENDIX C 


INTERACTION ENERGIES 


n 


l 


Since the essence of most long-range molecular approximations is the neglect of 
overlap between unlike orbitals (differently oriented p-states, for example) on colliding 
atoms, so that each a- or 7 r-bond or antibond exists more or less independently of the 
others, only the interactions which arise from a given pair of similar atomic functions 
are considered in this appendix. All the pertinent behavior is deduced from a study of 
the one-, two-, three-, and four-electron occupations associated respectively with H 2 , 
H 2 , He 2 , and He 2 ; in particular, theoretical relations between the different interaction 
energies are derived and compared with several techniques currently employed in the 
literature. 

Integrals are expressed in the Dirac "bra-ket" notation defined by 


<ip\p' • • • • *> s y [0(1)0’( 2 ) • • -]*FjV" (l)^” T (2) • • •]dT 1 dT 2 . . . 

(Cl) 

for the operator F and electrons 1, 2, • • • . 

The Coulomb and exchange energies Qj and Jj for an i-electron system are 
defined in terms of the attractive and repulsive interaction potentials and 4> r in 

the following manner: 


and 


so that 


Qi = 2 ^0a + 0r^ 

(C2) 

Ji = 2 " 1 (0a " 0r) 

(C3) 

0a, r = Qi ± J i 

(C4) 


A simple energy calculation involving the molecular orbitals a + b and a - b, 
where a and b are 1 s-atomic functions centered on the respective nuclei a and b, 
yields for the hydrogen molecule-ion H 2 the well-known expression 


0a = 


<a|H|a> + <a|H|b5 
1 + S 


-E 0 = 


<a]H|a>(l + S) + <a|H|b> - S<a|Hja> 


1 + S 


- E f 


= <a H a> - E 


o 


<a|Hlb> - S<a|H|a> _ _ 

+ Qi + Ji 

1 + S 


(C5a) 
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and similarly, 


, , <a|Hlb> - S<a|H|a> _ 

<p r - <a|H|a> - E 0 — Ql - 


1 - S 


(C5b) 


Other symbols represent the full Hamiltonian operator H, the ground-state energy E 0 
of atomic hydrogen, and the overlap integral S = <a|b>. 

Equations (C2) and (C3) then give 


Qj — <a H a> — Eq — SJ^ ~ e 2 R ^ — e^^ 


.-l| 


a> ~ 0 


(C6) 


and 


J, = (l - S2)' ^<a|H|b> - S<ajH|a>) * e 2 S<a 


r ' 1 

r b 


a> - e^<a 


r^ 1 b> 


(C7) 


terms of order S 2 being neglected to the right of the approximate equalities in order to 
depict large internuclear separations R. The potential -e2rg 4 refers to the interaction 
between the electron and nucleus b. 

An exactly analogous manipulation for H 2 , but this time with the valence-bond wave 
functions a(l) b(2) ± b(l) a(2) instead of molecular orbitals since the configuration inter- 
action of appendix A destroys the latter concept at large R, yields 


Q 2 = <ab|H|ab> - 2E 0 - S 2 J 2 ~ e 2 <ab r^ ab> - 2e 2 <a|rj" 1 |a> + e 2 R _1 « 0 (C8) 


and 


J2 


= (1 - s 4 ) 


(<aa|H|bb> - S' 


<ab|H|ab>) ~ 


e 2 <aa 


-1 

7 12 


bb> 


3 2 S 2 <ab 


-1 
7 12 


ab> + 2SJ 1 (C9) 


It is noted that whereas Jj is of order S, J 2 approaches zero as S 2 . This 
particular feature of long-range forces has been overlooked often in previous applications 
but is essential for subsequent comparisons. Such results, moreover, must not be con- 
fused with the statement (ref. 10, p. 362) that a two-electron bond is about twice as strong 
as the one-electron kind, for that conclusion is reached only by considering the near- 
equilibrium separations. At comparatively small values of R, the configuration inter- 
action does not so disrupt the concept of a doubly occupied molecular orbital in the ground 
state of H 2 ; the triplet repulsive state, on the other hand, is still described by 
a(l) b(2) - b(l) a(2) because 


18 



APPENDIX C 


\ 
f 

j a(l) + b(l) a(l) - b(l) a(l) b(l) 

2 a( 2 ) + b( 2 ) a( 2 ) - b( 2 ) a( 2 ) b( 2 ) 

The electron wave function for He 2 may be written 

[a(l) + b(l)]o<l) [a(l) - b(l[|a(l) [a(l) ± b(lj/3(l) 

|a(2) + b(2)]o<2) [a(2) - b(2)]o<2) [a(2) ± b(2)]/3(2) 

[a(3) + b(3)]c*3) [a(3) - b(3)]o<3) [a(3) ± b(3)]/3(3) 

a(l)o<l) b(l) a(l) [a(l) ±b(l)]/ 3 (l) 

= -2 a( 2 ) o< 2 ) b( 2 ) o< 2 ) [k( 2 ) ± b( 2 )]/ 3 ( 2 ) (Cll) 

a(3) o<3) b(3) a( 3) gi(3) ± b(3^/3(3) 

where a and j 8 refer to the two spin orientations. Thus, as was shown first by 
Linnett (ref. 12), only one molecular orbital exists, the remaining two electrons occupying 
atomic functions (AO). 

Since the AO electrons engage only in two-body exchange, which is of order S 2 , it 
seems reasonable to assume for present purposes that their sole function is to screen 
perfectly (by one electron unit each) the nuclei from the MO electron. Accordingly, one 
should take the a and b of the latter to be hydrogen 1 s, drop all reference to the AO 
electrons except as they provide proton-like cores for the MO, and write J 3 ~ Jj. This 
procedure is supported also by the fact that the MO electron is likely to be found at dis- 
tances farther from a nucleus than are those of atomic helium; consequently, the screening 
by the AO electrons is more effective in the molecule-ion. 

A similar study of the electron wave function for He 2 reveals that no molecular 
orbital is formed, even at small internuclear separations. If the atomic functions a 
and b are defined so that 

(t 1 + T 2 + e 2 r jg - 2e 2 r“* " 2 e 2 r 2a) a ^ = E ( He ) a ^ 2 ) ( C12 ) 

where Tj is the kinetic energy operator of electron i, there results 

$4 = Q4 “ J4 

Q 4 = <aabb|Hjaabb> - 2E(He) = 4e 2 <ab r~ 2 |ab> - 8 e 2 <a r^|a> + 4e 2 R“l ~ 0 (C13) 
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and 


J4 ~ 2^<aaab|H|abbb> - S 2 <aabb | Hj aabb>) = 4e 2 S<aa r ^ ab> + 2e 2 <aa 


r I2 bb> 


6e 2 S 2 <ab 


r j 2 ab> + 4e 2 s|s<a 2r^* a> - <a 


2r ; 


b> 


(C14) 


Contributions of order are neglected. 

The substitution into equation (C 14) of the approximate relation 


2<aa 


-1 

r 12 


ab> ~ S <aa 


-1 

r 12 


aa> + <ab 


-1 

r 12 


ab>, 


(C15) 


the foundation for which lies in the reasoning that a(i) b(i) corresponds to electron i 
dividing its time equally between a and b, gives 


J4 ~ 2e 2 <aa 


-1 

r 12 


bb> + 2e 2 S^aa 


-II 

' 12 | 


aa> - 2<ab 


-1 

'12 


ab>) + 4e 2 S|S<a 


2r ; 


-1 


a> - <a 


2r 


-1 


b> 


(C16) 


As illustrated next, this form is a step closer to the application of screening arguments 
similar to those for He^. 


Because of the saturation of atomic orbitals in the present example, the neglect of 
terms of order in equation (C14) refers strictly to pairwise exchanges of electrons 
with parallel spins. If each such pair is assumed to simultaneously create an antibond 
and perfectly screen the nuclei as seen by the other pair, only the interactions between 

electrons with like spins will survive the setting of a and b in equation (Cl 6) equal to 

9-1 2-1 

hydrogen 1 s functions. In addition, 2e^r^ must be replaced by e r^ to yield 


J4 ® 2e 2 <aa 


.-1 

12 


bb> - 2e 2 S 2 <ab 


,-l 

12 


ab> + 4SJj = 2J2 


(C17) 


A second way of reaching this same conclusion follows from the approximation of 
perfect pairing (APP) technique of Coulson (ref. 1) for the case of random spins. Since 
each electron of atom a pairs with each electron of atom b with one-fourth probability 
of antiparallel spins and three-fourths probability of parallel spins, there results 


04 = 4 




= -2J 2 


(C18) 
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In summary, one can write 


0i = Qi ± 


(i = 1, 2, or 3) 


$4 ~ Q4 “ J4 


Qi “0 


J3 ~ Jj = e 2 (S<a 


-1 

1 

-1 

r b 

a> - <aj 

r b 


b> 


and 


J 4 ~ 2J 2 = 2|e 2 <aajr j2|bb> - e 2 S 2 <ab 


,-l 

12 


ab> + 2SJ X I 


(C19) 


(C20) 


at large internuclear separations. As mentioned previously, a significant feature is the 
fact that the first pair J3 and Jj of equation (C19) bears no simple relation to the sec- 
ond pair J4 and J2 of equation (C20); it is only when equilibrium results are extrapo- 
lated incorrectly to large distances (without accounting for the effects of configuration 
interaction upon MO theory) that <p 2 85 2$j is obtained. Whereas the two-electron bond 
is stronger than the one-electron variety at short range, the roles are reversed further 
out. 

The previous statements regarding the existence of only one molecular orbital in 
the case of He^ and the failure of any to form for He2 require additional explanation. 
An illustration of the pertinent facts is obtained most conveniently by means of the alter- 
nate forms of the two-electron determinant of equation (CIO), according to which there 
are three methods of description as follows: 

(a) Molecular orbital description: Apart from the net electron-electron interaction 
(assumed to be negligible for this discussion) in J2, the assignment from equations (C4) 
and (C7) and the first form of equation (C6) of energy Jj(l - S) to the (a + b) electron 
and energy -Jj(l + S) to the (a - b) electron gives a total of -2SJ 4 ® -J2 for the two- 
electron antibond. Such additivity of independent one-electron energies is, of course, 
neither required nor implied by equation (CIO), but is imposed in order to take full advan- 
tage of the molecular orbital notion. 

(b) Valence bond description: This is the method of equations (C8) and (C9), which 
yields directly the antibond energy -J2. 

(c) Alternate molecular orbital description: The arbitrary assignment of energy 
J2/2 to the (a + b) electron and energy -3J2/2 to the (a - b) electron yields the correct 
antibond energy -J2 without resort to the -SJj portion of Qj used in description (a). 
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No allowance is made for the fact that J2/2 and Jj differ by at least the factor S. 
This technique is employed by Vanderslice, Mason, and Maisch in their work on nitric 
oxide (ref. 13) and oxygen (ref. 14). 

Of the three procedures, description (b) is preferred because there is no tendency 
to ignore electron-electron interactions by pursuing too literally, as in (a), the concept of 
independent one-electron energies. Consequently, there is less danger in the assertion 
that molecular orbitals do not exist in such examples. Description (c), on the other hand, 
although equivalent to (a) in its results for <p r (HjjJ and 0(He 2 j — visualize 

<#>( He 2 ) = -J 4 = 2(j 2 /2) + 2(-3J 2 /2) = -2J 2 

fails completely when the number of electrons is odd. More precisely, when one singly 
occupied MO is the essential feature, as in 0 a (H e 2 )> and 4> r (He'£j, 

description (c) gives the incorrect answers J 2 /2, -3J2/2, -J2/2, and -5J2/2, respec- 
tively. Each of these results disregards the factor S”* and leads to serious diffi- 
culties, among which is the loss of the attractive character of $ a ^He^j. In addition, 
asymmetries are introduced between the bonds and antibonds of the same molecule. 

Finally, should the electron-electron interactions of equation (C20) be comparable 
to 2SJj, the literal interpretation of description (a) will be insufficient. A typical 
numerical calculation (see ref. 10, pp. 329, 342, and 343, for the appropriate integrals) at 
R = 4 atomic units shows that the neglected terms in (a) constitute about 40 percent of 
those retained; accordingly, Jj must be considered in at least the complexity suggested 
by equations (C19) and (C20). More elegant treatments, such as the generalization by 
Kolos and Wolniewicz (ref. 15) of the James and Coolidge wave function to large distances 
(keeping the dependence upon rj 2 ^, are required for absolute accuracy of the detailed 
characteristics (including the Van der Waals forces) of the potential tails. 

These discussions, and especially the overlap distinctions between one- or three- 
and two- or four-electron configurations, are particularly significant for molecules in 
which two or more different types occur simultaneously. Previous calculations on the 
long-range interactions of nitrogen and oxygen (ref. 13) and two oxygen atoms (ref. 14), 
for example, must be reevaluated because of the relations employed between the two- and 
three-electron bonds and antibonds existing together in certain states. 
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